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Spanning avalanches in the 3D Gaussian Random Field Ising Model (3D-GRFIM) with metastable 
dynamics at T = have been studied. Statistical analysis of the field values for which avalanches 
occur has enabled a Finite-Size Scaling (FSS) study of the avalanche density to be performed. 
Furthermore, direct measurement of the geometrical properties of the avalanches has confirmed an 
earlier hypothesis that several kinds of spanning avalanches with two different fractal dimensions co- 
exist at the critical point. We finally compare the phase diagram of the 3D-GRFIM with metastable 
dynamics with the same model in equilibrium at T = 0. 
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I. INTRODUCTION 

Avalanche behavior has been found in many first 
order phase transitions at low temperatures. When the 
external parameters are slowly driven, the transition, in- 
stead of occurring as a sharp change of the system prop- 
erties at a certain point of the phase diagram, splits into a 
series of discontinuous jumps that link metastable states. 
This gives rise to hysteresis loops, whose branches consist 
of a sequence of random steps. In many cases the sizes of 
such steps range from microscopic to macroscopic scales, 
distributed according to a power-law. The phenomenon 
has been found to be associated with magnetic transi- 
tions 0, 0, 3 J capillary condensation 0, U, martensitic 
transformations 0, and others @. A crucial ingredient 
in order to observe such avalanches is that thermal fluctu- 
ations are very small compared with the energy barriers 
that separate transformed and untransformed domains. 
For this reason, such first-order phase transitions have 
been called "athermal" 01 or "fluctuationless" [TollTTI |. 

Within this context the 3D-GRFIM with metastable 
dynamics is a prototype model: the complexity of what 
we call "disorder" in a real system is simplified into a 
series of quenched random fields, Gaussian distributed 
with zero mean and standard deviation cr, that act on 
every spin of a 3D Ising model. In addition, one assumes 
that temperature is zero (T = 0) and provides the model 
with a particular metastable dynamics in order to study 
the evolution of the magnetization m when the external 
field H is swept. The details of the dynamics were in- 
troduced by Sethna and coworkers a decade ago 0, 0| . 
The basic assumption is that the driving field rate is slow 
enough so that system relaxation can be considered in- 
stantaneous (adiabatic driving). Such relaxations are the 
so-called magnetization avalanches. 

After the introduction of the model several works 



El 111 111 mill in described the basic associated 
phenomenology: (a) the existence of a disorder-induced 
critical point at <t = Uc associated with the change from 
a continuous to a discontinuous hysteresis loop and (b) 
the fact that within a large region around the critical 
point the distribution of avalanche sizes D[s) exhibits 
almost power-law behavior. Nevertheless, several ques- 
tions remain unsolved, mostly related to the properties of 
the spanning avalanches which are responsible for the ob- 
served macroscopic discontinuities in the H-m hysteresis 
loop and in the a-H phase diagram. 

More recently, a FSS analysis of the number of 
avalanches and avalanche size distribution |20l | has re- 
vealed that the scenario is quite complex. The study 
was restricted to the statistical analysis of the full set 
of avalanches recorded in a half loop, irrespective of the 
field values H where such avalanches occur. (The ob- 
tained statistical distributions are often called integrated 
distributions). A detailed study as a function of H has 
not been done before. 

Let us summarize here the main results presented in 
Ref. I23, in order to introduce the notation. Accord- 
ing to the dependence on the system size L of the av- 
erage number A^q (cr, L) and size distribution Da{s; a, L), 
avalanches can be classified into several categories, (the 
subindex a stands for the different categories) as pre- 
sented in Table m The scaling behavior is written as a 
function of the variable u which measures the distance 
to the critical value of the disorder CTc = 2.21 ± 0.02. Its 
precise definition is given as a second-order expansion: 
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with A = —0.2. This expression was found to be the 
best choice for the collapse of the scaling plots. The 
values of the different critical exponents are summarized 
in Table ITTl together with the new exponents that will be 
computed in the present work. 

Classification of the avalanches starts by checking 
whether the avalanches span the system in 1,2 or 3 spatial 
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avalanche type a average number size distribution 

non-spanning ns N„s{a,L) [N„s = N„sc + N„so] D„s{s;cr, L) [N„sD„s = N„scD„sc + N„soD„so] 

critical non-spanning nsc Nnsc{(7,L) = L""--" Nnsc(uL^/'') Dnsc{s;a,L) = L"'""''^f D3c{sL~'^f ,uL^^'") 

non-critical non-spanning nsO N„so{o', L) = L"^ N„so{o') 

ID-spanning 1 Nija, L) ^ L" Ni(uL'^^") Di{s;o, L) = L-'^f Di{sL-'^i ,uL^'") ~ 

2D-spanning 2 Afaja, L) = L^jVai^L^''") D2 {s;o,L) = L-"^! Di jsL-'^i ,uL^'") 

SD-spanning 3 A^aC^.i) [iVa = A^Sc + iVs-] D3(s;cr,L) [iVsDa = A^ScI^sc + A^a-Da-] 

critical 3D-spanning 3c Nz^{a,L) = L" N3^{uL^''') Dsds-a, L) = L-''f DacisL-'^f ,uL^^'') 

subcritical 3D-spanning 3- Ng-ja, L) ^ N-j-juL^^") D3-{s;a,L) = L'"^^- D^-jsL^'^^- ,uL^/'') 



TABLE I: Classification of the different types of avalanches in the 3D-RFIM according to their geometrical properties and their 
finite-size scaling behavior. Non-spanning avalanches and 3D-spanning avalanches exhibit mixed scaling behavior (indicated in 
square brackets). They can be separated into two subcategories. 
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2.78 ±0.05 


2^, this work 
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2(^, this work 
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0.15 ±0.08 






0.024 ±0.012 


20 


1/m 


1.5 ±0.1 


this work 



parameter using (3c = j^(3 — 9 — df) = 0.15 ± 0.08. 
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TABLE IL Summary of the values of the critical exponents of 
the 3D-GRFIM with metastable dynamics, obtained in Ref. 
l2fll and in the present work, as indicated in the last column. 



directions. |33l| These three kinds of spanning avalanches 
are indicated by a = 1, 2, 3 respectively. 

To obtain good scaling collapses of the 3D-spanning 
avalanches, an extra hypothesis was introduced: they can 
be separated into two subcategories, subcritical (a = 3—) 
and critical {a = 3c) that scale with different exponents. 
In particular, this assumption indirectly leads to the 
conclusion that they must have different fractal dimen- 
sions (is- and df. Moreover, non-spanning avalanches 
(a = ns) should also be separated into two subcategories: 
non critical (a = nsO) and critical {a ~ nsc), depend- 
ing on whether their number and size distribution scales 
with distance to the critical point u or not. 

Fig.^presents a summary pH of the scaling functions 
Na according to the results in Ref. From the be- 

havior of such scaling functions when uL^^" — > ±oo one 
can sketch out the scenario in the thermodynamic limit. 
Below <Tc one subcritical 3D-spanning avalanche exists, 
which is responsible for the observed discontinuity of the 
magnetization in the thermodynamic limit. Such a dis- 
continuity is the order parameter and vanishes when ap- 
proaching f7~ according to Am oc (<Tc - cr)'^3- ^ with an 
exponent /33_ = v{?, - dg^) = 0.024 ± 0.012. Never- 
theless, it is difficult to find this critical behavior from 
simulations of finite systems since the contributions from 
critical spanning avalanches (a = 1,2 and 3c) may lead to 
quite good (but incorrect) scaling collapses of the order 
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FIG. 1: Scaling functions corresponding to the number of 
ID-, 2D-, critical 3D- and subcritical 3D-spanning avalanches 
in the 3D-RFIM as found in Ref. The symbols show the 

overlap of the data corresponding to numerical simulations 
with many different system sizes ranging from L = 5 to L = 
48. 

A large number of non-critical non-spanning 
avalanches exist for the whole range of a. They 
cannot contribute to any observed macroscopic jump 
since their size is vanishingly small in the thermody- 
namic limit. At cr = (Tc the 6 categories of avalanches 
exist. On average, one finds 0.79 ± 0.02 subcritical 
3D-spanning avalanches and an infinite number of the 
five other types of avalanches. At the critical point (or 
close enough to it), the distribution of avalanche sizes 
is dominated by critical non-spanning avalanches and 
exhibits an approximate power-law behavior with an 
exponent Te// = r„sc + (3 — Onsc)/df = 2.00 ± 0.06. This 
power law behavior is restricted to the central part of 
the avalanche size distribution 1 <^ s <^ . 
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In the present paper we will concentrate on the 4 differ- 
ent kinds of spanning avalanches, extending the finite-size 
scaling analysis and focusing on (i) the study of the val- 
ues of the external field for which the spanning avalanches 
occur and (ii) the direct measurement of the geometrical 
properties of the avalanches. This will enable the critical 
exponent /i, related to the renormalization group (RG) 
flow along the field direction, to be found and have a 
direct test of the FSS hypothesis that leads to the sepa- 
ration of the two kinds of SD-spanning avalanches (crit- 
ical and subcritical) with different fractal dimensions: 
df = 2.78 ± 0.05 and dg. = 2.98 ± 0.02. 

In section ^ we summarize the 3D-GRFIM and the 
details of our numerical simulations. In section ITTll raw 
numerical results are presented. In section IIVI we dis- 
cuss the main FSS hypothesis, as an extension of those 
presented previously. These hypotheses are checked in 
section El In section |Vll we determine the geometrical 
properties of the avalanches. In section IVIII we discuss 
the consequences of the present study and, finally, in sec- 
tion ^^H] we summarize and conclude our findings. 



II. MODEL 

The 3D-GRFIM is defined on a cubic lattice of size 
L X L X L with periodic boundary conditions. On each 
lattice site (i = 1,...,L^) there is a spin variable Si 
taking values ±1. The Hamiltonian is: 



(2) 



i=l 



where the first sum extends over all nearest-neighbor 
(n.n) pairs, H is the external applied field and hi are 
quenched random fields, which arc independent and are 
distributed according to a Gaussian probability density 
with zero mean and standard deviation a. 

The equilibrium ground-state (T = 0) of this Hamilto- 
nian has been recently studied |22l | . In this work we focus 
on the metastable version of the 3D-GRFIM proposed for 
the analysis of the behavior at T = when the system is 
driven by the external field H. For H = +oo the state 
of the system which minimizes Ti. is the state with max- 

imum magnetization m = ^i^i Si/ L = 1. When the 
external field H is decreased, the system evolves follow- 
ing local relaxation dynamics. The spins flip according 
to the sign of the local field: 



H 



(3) 



where the sum extends over the 6 nearest-neighboring 
spins of Si. Avalanches occur when a spin flip changes 
the sign of the local flcld of some of the neighbors. This 
may start a sequence of spin flips which occur at a fixed 
value of the external field H , until a new stable situation 



is reached. H is then decreased again. The size of the 
avalanche s corresponds to the number of spins flipped 
until a new stable situation is reached. Note that the 
corresponding magnetization change is Am = 2s/ L^. 

The numerical algorithm we have used is the so-called 
brute force algorithm which propagates one avalanche at 
a time We have studied system sizes ranging from 
L = 5(L3 = 125) to L = 180(L3 = 5832000). The mea- 
sured properties are always averaged over a large num- 
ber of realizations of the random field configuration for 
each value of ct, which ranges between more than 10^ for 
L < 80 to 300 for L = 180. 

We have recorded the sequence of avalanche sizes dur- 
ing half a hysteresis loop, i.e. decreasing H from +oo 
to —OQ. We have determined not only the size s of 
each individual avalanche, but also the field H at which 
each avalanche occurs. The avalanches have been clas- 
sified as non-spanning, ID-spanning, 2D-spanning and 
3D-spanning as explained in Ref . (2Q] . 

By performing statistics, we obtain the average density 
of avalanches for each type occurring within an interval 
{H, H + dH) [111 . We will call this quantity the number 
density na{H] a, L). It satisfies: 



dH na{H;a,L) ^ Nai(T,L) 



(4) 



where Na (g, L) are the average number of avalanches of 
each kind [2fl| defined in Table We have also mea- 
sured the bivariate size distribution T^ais, H\ a, L). This 
probability density is normalized so that: 



/ dHVa{s,H;a,L) = 1 

1 ^ —oo 



(5) 



Note that by projecting Pq, we can obtain the two 
marginal distributions: 



,ff;a,L) 



na{H;a,L) 
N^ia,L) 



(6) 



and 



/oo 
dHV^{s,H-a,L)^D^{s-a,L) (7) 

which represent the probability density of finding an 
avalanche of type a within [H, H + dH) and the proba- 
bility (integrated distribution) that an avalanche of type 
a has a size s (after the half loop) respectively. 

In the following sections bivariate distributions will be 
presented as point clouds and the marginal distributions 
as histograms. Point clouds provide qualitative under- 
standing of the distributions. Quantitative analysis is 
much better performed from the marginal distributions 
and their moments. The average field where the different 
types of avalanches occur will be particularly interesting: 



dH H 



na{H; a, L) 



(8) 
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and its standard deviation: 



,1/2 



(9) 



In some simulations we have also performed a sand 
box counting analysis j23 | of each individual spanning 
avalanche in order to have a direct determination of their 
fractal dimension. The method consists of considering 
boxes of linear size £ (ranging from 1x1x1,3x3x3, 
5 X 5 X 5, ... up to (L — 1) X (L- 1) X (L- 1)) centred near 
the spin that triggered each spanning avalanche. We de- 
termine the number of sites that belong to the avalanche 
for each box. By averaging over all the avalanches of the 
same kind (occurring during a half loop) and over many 
disorder realizations, we obtain the mass Ma{i]cr^L) as 
a function of the box length £. 



III. NUMERICAL RESULTS 

Fig. 121 shows a point cloud corresponding to 
2?3(s, H; a, L) for L = 48 and different values of a. As can 
be seen, below (Tc 3D-spanning avalanches are large, close 
to s = and exhibit a certain spread around a value 
of H which shifts upwards when CTc is approached from 
below. At cr = (Tc the avalanches concentrate around 
a critical field value He — —1.42 and span almost all 
possible sizes. Above Cc the few existing 3D-spanning 
avalanches are small. In fact, they are smaller than the 
mean size of a percolating cluster in a 3D cubic lattice 
(s ^ 0.311L3). 
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FIG. 2: Point clouds corresponding to SD-spanning 
avalanches for a system with size L — 48. Each point cor- 
responds to an individual avalanche and indicates the size s 
and the field where it occurred H. Clouds corresponding to 
different values of a are indicated by different symbols. Data 
have been averaged over many realizations of disorder. The 
dashed line indicates the maximum value s — L^. 

Fig. 121 shows the point clouds corresponding to non- 
spanning, ID-, 2D-, and 3D-spanning avalanches at a — 



Gc- Non-spanning avalanches are suppressed from the 
main plot, but are shown in the inset on a log-linear scale. 
As expected, spanning avalanches concentrate around 
He- 3D-spanning avalanches exhibit a larger size than 
2D-spanning avalanches and the latter show a size larger 
than ID-spanning avalanches. Moreover, one can see that 
3D-spanning avalanches are distributed in a double cloud. 
As will be seen, these two clouds correspond to the two 
types of 3D-spanning avalanches predicted in |20j . 
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FIG. 3: Point clouds corresponding to ID-, 2D-, and SD- 
spanning avalanches at cr = CTc for L = 48. The inset shows, 
on a log-linear scale, the same data together with the cloud 
corresponding to non-spanning avalanches. 

The dependence on the system size L is illustrated in 
Fig. 21 The point clouds represent the distribution of 
spanning avalanches at ct = (Tc for increasing values of L. 
The three plots correspond to (a) ID-, (b) 2D-, and (c) 
3D-spanning avalanches. The dashed line corresponds to 
our estimated value of He — —1.425. Note that all kinds 
of spanning avalanches tend to concentrate around such 
a field value for increasing system sizes. The shape of the 
cloud corresponding to 3D-spanning avalanches remains 
assymetric, while L grows. 

To obtain a quantitative measure of H,., we have com- 
puted the field averages {H)2 and (ff)3. Their be- 
havior as a function of L at cr = tJc is plotted in Fig.[SJa). 
In the next sections a FSS analysis will be formally pro- 
posed. Nevertheless, from the behavior in Fig. E^a) one 
can already guess the following scaling hypothesis: 



(10) 



where /i will be the exponent governing the divergence of 
the correlation length when the field H approaches He at 
a = Uc- A first check of this hypothesis is performed in 
Fig-E^b) by a 3-parameter {Hc^ and Ca) least-squares 
fit to the three sets of data. We have obtained the same 
values of 1/^ = 1.5 ± 0.1 and H^ = -1.425 ± 0.010 for 
{H)i, {H)2 and {H)^. Results of the fits are indicated 
by dotted lines in Figs.[S{a) andE^b). 
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FIG. 4: Point clouds corresponding to ID- (a), 2D- (b) and 
3D-spanning avalanches (c) at a = ac for increasing system 
sizes, as indicated by the legend. The same scale is used on 
the three plots. 
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IV. FINITE-SIZE SCALING HYPOTHESIS 

A. Scaling of distributions 

To proceed with the analysis of the numerical data one 
must postulate the ad-hoc FSS hypothesis. As done in 
Ref. |23|, we follow standard RG arguments. The dis- 
tance to the critical point [ucHc) is measured with two 
scaling variables u{a,H) and v{a,H)^ which depend on 
the externally tunable parameters a and H . After a 
renormalization step that transforms the system size L 
as: 



Ly, = b-^L. 



(11) 



u,v, and the size s of an avalanche of type a will behave 
as: 



Ub 



U Vb 



V Sb 



(12) 



where da are the fractal dimensions of the avalanches 
and the exponents v and fi control the divergence of the 
correlation length when approaching the critical point 
along the two directions u and v. Consequently, in order 
to formulate the FSS hypothesis, we will consider the 



FIG. 5: (a) Average fields {H)a{(Jc, L). The dashed lines indi- 
cates He ~ —1.425 and the dotted lines correspond to power- 
law fits, (b) The same data represented on log-log scales. The 
lines indicates the power-law behavior. Note that the 3 kinds 
of spanning avalanches exhibit the same exponent. 



following three main RG invariants: 



vL 



(13) 



(14) 



(15) 



Therefore, close enough to the critical point the bivariate 
distributions T^ccis, H] a, L) will behave as: 



1/m 



(16) 

where is df ~ 2.78 for the ID-, 2D- and critical 3D- 
spanning avalanches and d^- ~ 2.98 for the subcritical 
3D-spanning avalanches, as found in Ref. [20j . 
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The number densities na{H; a, L) will behave as: 

n^{H; cr, L) = L^^+i/^ria (uL^/", vL^/^'^ , (17) 

where 0a is 0c = 0.1 for the ID-, 2D- and critical 3D- 
spanning avalanches and 9a = for the subcritical 3D- 
spanning avalanches. The FSS hypothesis l|16|l and H17() 
extend the hypothesis presented in Tabled including the 
dependence on the external field H. 

B. Detailed scaling variables 

The proper dependence of the scaling variables u and 
V on the tunable model parameters a and H is unknown, 
but it should be analytic 0il23- Keeping this in mind, 
we consider the following second-order expansion around 
a — <7c and H = He'- 

u = ui+Aul + A'vi+A"uivi+A"'vl (18) 
V = vi+Bvl + B'ui + B"uivi+B"'ul, (19) 

where ui = {a — (Jc)/<^c and vi = {H — He) /He are the 
first-order scaling variables. These variables represent 
the simplest choice to measure the distance to the critical 
point. Now we want to check which corrections to the 
first-order scaling variables arc really important for large 
system sizes. Consider a function F{H; a, L) that can be 
written as 

F{H; (T, L) = L'^FiuL^/", vL^'^), (20) 

when the correct scaling variables are considered. In the 
case in which we try to obtain scaling collapses with two 
incorrect scaling variables u' and w', an explicit depen- 
dence on L may appear, which makes the collapse L~°'F 



impossible for different system sizes. In such a situation: 

F(iJ; a, L) = L''G{u' L^''' , ^'L^/^, L), (21) 

where we have introduced a "non" -scaling function G 
that depends explicitly on L. By construction, G coin- 
cides with F{uL^^'^ , wL^/^) when u' = u and v' = v. The 
effective scaling variables u' and v' can be expanded up 
to second order: 

u' = ui + Au\ + A'vi + A"uivi + A"'vl (22) 
v' = vi+Bvl + B'ui+B"uivi+B"'ul. (23) 

Up to second order, the difference between 
G{u'L^/'',v'L^/^',L) and FiuL^'" ,vL^/ ^') is 

Giu'L^'", i/L^/^, L) - FiuL^'", vL^/'') = 
L^'"" Fi{uL^'\vL^'^'){u' ^ u) 
+ L^/^'F2{uL^''',vL^/^'){v' ^v) 

+ Li/'^+i/^ A2(wii/", vL^'^'){u' - u){v' - w),p4) 
where the subindex 1 and 2 in the scaling functions F in- 
dicate the partial derivatives with respect to uL^^" and 
vL^^f^, respectively. Introducing the expansions 118|l . 
(Unj, (1211), and |(221) in Eq. defining an appropri- 

ate set of eleven functions /^(uiL^/'^, wiL^/^) from the 
derivatives of F, and re-ordering the terms in (|24|) we 
can write: 



G(u'Li/^ u'Li/^, L) - F(uLi/^ vL'/^) = {A' - A')L'/''-'/^ (^) A + (B' - B')L'/^-'^-' (^) ^ 
+ [A- A)L-^/^ (^) ' h + {B- B)L-'/^ (^) ' U 
+ (A" - A")L-y^\^) (^) .h + {B" - B'OL-V. A 
+ {A!" - A"')L^'''~^'^ + - B"')L''^-^''' (^)' 

+ (^' - A'fL^/--^/^^ ( v) ' + " 5')'^'/^"'/" (^) ' ho 

+ (^'-A')(B'-S')(^)(^)/ii, (25) 



where the dependence of the functions fi on ui L^/" and 
ViL^^^ is not written for simplicity. 



Given that l/i/ 0.8 (Tablc|lTI) and 1/^ = 1.5 (Fig.Ell, 

and noting that only the terms multiplying positive pow- 
ers of L will be important in the thermodynamic limit. 
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we conclude that only terms in which B' — B' appears 
represent an exphcit dependence of G{u' L^^" ,v' L^^ ^ , L) 
on L. Therefore, only the B' term is relevant in the 
thermodynamic limit and thus must be considered in the 
expansions H22fl and H23() . The remaining coefficients may 
be neglected. In particular, a term proportional to u\ is 
not necessary if we consider large system sizes. Such an 
irrelevance was in fact observed in Fig. 8 of Ref. [20l |. 
Nevertheless, we will retain the quadratic correction Au\ 
in order to compare appropriately with previous results 
|20l |. In summary, we will use the following approxima- 
tions for the scaling variables: 

u = ui+ Aul, (26) 
V ^ VI +B'ui. (27) 

The correction (|27|l associated with the distance to He 
was first introduced in Ref. [2^ where the parameter 
analogous to B' was called the "tilting" constant. We 
should mention that the authors demonstrate the im- 
portance of such a correction with arguments that are 
slightly different those proposed here. 

C. Scaling of the field averages and standard 
deviations 

One is now ready to deduce the scaling behavior cor- 
responding to the field averages defined in Eq. © and 
to the standard deviations defined in Eq. I^. On the 
one hand, multiplying the marginal distribution Ua/Na 
by iJ, integrating over the full H range, and using the 
relation H27|) . we find 

(1 - B'ui) - (if)„(a, L) = L-^/^'K {^L^'") ■ (28) 

It is useful to define an "effective" disorder-dependent 
critical field H*{(t) as: 

H^ia) = Hcil- B'ui). (29) 

For (T = tJc we recover the scaling hypothesis proposed 
in Eq. H1U|I with Cq = /iq(0). In this case we obtain an 
estimate of l//.t that is unaffected by the tilting constant 
B'. 

On the other hand, by performing similar calculations, 
it is easy to write cr^ {a, L) as: 

a^{a,L)=L-''^a^{uL''''). (30) 

Notice that this scaling expression is also unaffected by 
the tilting constant. 

D. Separation of the two kinds of 3D-spanning 
avalanches 

From the FSS analysis of the integrated distributions, 
it was suggested in |2y that two kinds of 3D-spanning 



avalanches exist with different fractal dimensions. This 
assumption allowed for excellent collapses of the scaling 
plots. Nevertheless, the separation of the scaling func- 
tions corresponding to the two kinds of avalanches was 
possible by using a double FSS technique involving the 
collapse of data corresponding to three or more different 
system sizes. The propagation of the statistical errors 
within such complicated computations rendered large er- 
ror bars in the scaling functions and exponents. 

Given the two different fractal dimensions, it would 
be desirable to be able to perform a direct classification 
of the 3— and 3c avalanches during simulations. Nev- 
ertheless, this desirable idea is not possible since, in a 
finite system, a good determination of a fractal dimen- 
sion is only possible after performing statistics of many 
avalanches of the same kind. 

In this work we propose two separation methods 
that, although being approximate (a small fraction of 
avalanches are not well classified) give enough bias to the 
statistical analysis to allow for determination of the dif- 
ferent properties of subcritical and critical 3D-spanning 
avalanches. 

The idea behind the methods is that for a finite sys- 
tem that is below CTc, one basically finds one subcriti- 
cal 3D-spanning avalanche. The other types of spanning 
avalanches may occur only close to Uc- Moreover, given 
their different fractal dimension, we expect subcritical 
3D-spanning avalanches to be larger. Thus, we propose 
the following two methods, which will be applied only 
below (Tc- 

Method 1: The larger 3D-spanning avalanche in a 
half loop is classified as subcritical. The other 3D- 
spanning avalanches will be considered critical 3D- 
spanning avalanches. (We have checked that the larger 
3D-spanning avalanche is also the last 3D-spanning 
avalanche found when decreasing the field from H — -\-oo 
to H ~ —CO in almost all the studied cases.) 

Method 2: We classify a 3D-spanning avalanche as sub- 
critical only when no other spanning avalanches occur 
during the half loop. If other spanning avalanches occur, 
we classify them as critical SD-spanning. The idea behind 
this method (which we will discuss in section IVIip is the 
conjecture that the subcritical 3D-spanning avalanche, 
close to, but below (Jc fills a large fraction of the system 
and does not allow other spanning avalanches to exist. 

Table IIIII shows how the two methods classify a cer- 
tain 3D-spanning avalanche depending on whether the 
other spanning avalanches found in the half loop are ID- 
, 2D- or 3D-spanning. In this latter case the fact that 
the other 3D-spanning avalanche(s) found are smaller or 
larger than the avalanche being classified must be taken 
into account. The two methods only differ in the case in 
which the 3D-spanning avalanche being classified is the 
largest, but other spanning avalanches (either ID, 2D or 
smaller 3D) exist in the loop. 

Fig. shows an example of separation of 3D-spanning 
avalanches into subcritical and critical, using the two 
methods. It corresponds to i = 48 and ac = 2.21. One 
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ID, 2D or smaller 
3D-spanning 
avalanche exist 


larger 3D- 
spanning 
avalanche 
exist 


Method 1 


Method 2 


no 


no 


3- 


3- 


yes 


no 


3- 


3c 


no 


yes 


3c 


3c 


yes 


yes 


3c 


3c 



TABLE III: Summary of the classification of a given 3D- 
spanning avalanche according to the two methods proposed 
in the text. 



can appreciate that the original double-shaped cloud is 
separated into two. The cloud corresponding to critical 
SD-spanning avalanches is similar in shape to the clouds 
corresponding to ID- and 2D-spanning avalanches. Note 
that Method 2 classifies a certain number of large 
avalanches, occurring at very negative fields, as being 
critical that Method 1 classifies as being subcritical. 
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FIG. 6: Example of separation of 3D-spanmng avalanches 
into subcritical and critical, using methods 1 and 2 explained 
in the text. 

The two separation methods will be used throughout 
the rest of the text to separately analyze the data corre- 
sponding to subcritical 3D-spanning avalanches and crit- 
ical 3D-spanning avalanches. In some of the statistical 



analysis presented below we will consider only the SD- 
spanning avalanches which are equally classified by the 
two methods and discard those which are classified differ- 
ently from the analysis. Although this procedure reduces 
the size of the statistical sample, it ensures that we do 
not introduce any bias due to ill-classification of some of 
the avalanches. 



V. SCALING COLLAPSES 
A. Field Averages and standard deviations 

Fig. [71 presents the scaling collapses corresponding to 
the field averages for cr < CTc- Data is presented on log- 
log scales in order to analyze the power-law behavior for 
\u\L^^'^ oo. Fig. Efa) shows data corresponding to 
{H)i, {H)2, and {H)^^ whereas Fig. [7|[b) shows data 
corresponding to {H)^-. We remark that only the SD- 
spanning avalanches equally classified by Methods 1 and 
2 have been used for computing the averages correspond- 
ing to (i?)3- and (H) :ir-. By fixing l//i = 1.5 and tak- 
ing = 1.2 from Ref. [23| we get the best collapses for 
B' = 0.25 ± 0.10. We would like to emphasize that the 
four sets of data scale extremely well with the same val- 
ues of fj, and ly. 

The asymptotic behavior of /i3_ (dotted fine in Fig. 

Hb)) for large values of lulL^/" is S.3 (|u|Li/'') The 
exponent 1.8 equals i^/fj, within statistical error. This 
means that, in the thermodynamic limit. 



{H)3^{a) = H:ia)~3.3\u\''/^. 



(31) 



Fig. |S1 shows this behavior, which finishes at the crit- 
ical point [ac, He) because no subcritical SD-spanning 
avalanches exist above. The disorder-dependent criti- 
cal field H*{a) and the critical field He arc indicated 
by dashed and dotted lines, respectively. We have also 
plotted the numerical estimates of (i?)3- for different 
system sizes in order to show that Eq. (|S1|I is the lim- 
iting behavior for L — > oo. The asymptotic behavior of 
ha for the ID-, 2D- and critical SD-spanning avalanches 
is proportional to (|u|L^/'^) . This implies that in the 
thermodynamic limit 



{H)ai<7)=H:ia) 



(32) 



for the ID-, 2D- and critical SD-spanning avalanches. 

Similar finite-size scaling analyses can be done for the 
standard deviations of the marginal distributions ria/Na 
according to Eq. From the obtained collapses we 

deduce the following behavior for large values of \u\L^^'^ 
(cr < ac): d^_^{uL^/'') ~ {\u\L'^''')°-^ for the subcritical 
SD-spanning avalanches and (tU(mL-^/'') ^ (|m|L^/'^)°-^ 
for the ID-, 2D-, and critical SD-spanning avalanches. 
Similar behavior is observed for a > Uc- These results 
(see Eq. (|SO|l 'l indicate that the standard deviation of 
the marginal distribution Ua/Na corresponding to any 
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FIG. 7: Scaling plots corresponding to the difference between 
{H)a and the effective critical disorder Hc{cr) for (a) ID-, 2D- 
, and critical SD-spanning avalanches and (b) the subcritical 
SD-spanning avalanches. We present only data for a < ac. 
The dotted lines correspond in each case to the asymptotic 
behavior of —ha when ImIL^/" ^ oo. The scaling plots corre- 
sponding to the ID- and 2D-spanning avalanches have been 
displaced one decade upwards for clarity. 



type of spanning avalanche vanishes in the thermody- 
namic limit for any value of a. 



B. Number density 

The number density corresponding to the ID-spanning 
avalanches at the critical amount of disorder ni(H; L) 
is shown in Fig. Ela) as a function of the applied field 
for different system sizes. The number density shovifs a 
peak that increases and shifts for increasing L. Similar 
behavior is observed for n2{H; ac, L) (Fig.^Ja)). A FSS 
analysis is performed using the scaling assumption for the 
number densities (Eq. (|17|l '). The results of such an anal- 




FIG. 8: Representation of {H)-s-{cf) (Eq. JSTJ), Hl{a) 
(Eq. 129^ '). and He. Symbols correspond to the numerical 
estimate of (-ff)3~ from simulations for different L, as indi- 
cated by the legend in the previous figure. 



ysis are presented in Figs.^fb) and irUT bl for ni{H; ac, L) 
and n2{H;ac,L), respectively. To obtain these collapses 
we have used = 0.1 and = 1.5 as in the preceding 
sections. The scaling functions in Figs. Efb) and llOf b) 
are well approximated by Gaussian functions (indicated 
by continuous lines). When wL^/^ ±cx) both scaling 
functions go exponentially to zero. This behavior indi- 
cates that, in the thermodynamic limit for a = ac, ID- 
and 2D-spanning avalanches only exist at H = He- 

Fig.llllshows several cuts corresponding to the scaling 
surface ni{uL^/^ ,vV-/'^). From the collapses we obtain 
B' ~ 0.25 ± 0.10 in total agreement with previous esti- 
mates. From a qualitative point of view, the collapses 
indicate that the scaling surface shows a crest with am- 
plitude depending on uL}l" . More quantitatively, the 
scaling collapses for each cut can be well approximated 
by Gaussian functions, whose amplitude, peak position, 
and width depend on uL^^^ . Furthermore, the depen- 
dence on uL^/'^ of the fitted amplitudes also adjusts very 
well to a Gaussian function that follows the profile of the 
crest (continuous line on the back plane in Fig. lll|) . The 
dashed line on the bottom plane indicates the position 
of the crest {vL^/ t')i{uL^/'') = hi{uL^/'') / He, which has 
already been shown in Fig. [7Ja). 

All these considerations imply that, for any value of 
H, the scaling function hi decays exponentially when 
uL^/'^ — > ±oo. This indicates that, in the thermody- 
namic limit, irrespective of the value of H , ni is zero for 
a ^ ac- In contrast, when a = ac, n-i diverges at H = He 
and is zero for other values of the field. This scenario for 
ni is also applicable to n2. 

To obtain n^c and ns- we have used Method 2 of sep- 
aration described in Sec. IIVDI The results for a = ac 
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FIG. 9: (a) Number density of spanning avalanches in one 
dimension at the critical amount of disorder, (b) Scaling plot 
corresponding to the data in (a) according to Eq. 1171 1 with 
9 — 0.1. The continuous line in (b) shows a Gaussian fit. 
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FIG. 10: (a) Number density of spanning avalanches in two 
dimensions at the critical amount of disorder, (b) Scaling plot 
corresponding to the data in (a) according to Eq. 1171 1 with 
9 — 0.1. The continuous line in (b) shows a Gaussian fit. 



VI. DIRECT DETERMINATION OF THE 
FRACTAL DIMENSIONS 



are presented in Figs. [T^ and 1131 for ji^dH; ac, L) and 
n3-{H;ac,L), respeetively. We have also tried to sepa- 
rate the two kinds of avalanches by Method 1, but the 
collapses are not as good as those obtained with Method 
2. This result indicates that in the set of avalanches 
non-equally classified by the two methods, there are 
more critical SD-spanning avalanches than subcritical 
SD-spanning avalanches. 

As in the case of rii and 77,2 at cr = Cc, the behavior 
of the scaling functions in Figs. ll2f b'l and ll^f b'l indicate 
that both n^c and 713 _ diverge at H — He and are zero for 
fields different to He- The detailed study of the bivariate 
collapses corresponding to n^dH; a, L) and n^-^H; a, L) 
for a ^ ac is difficult and tedious. In particular, for a > 
CTc we do not expect the separation methods to work and 
for a < ac, the analysis would require a lot of statistics. 



In the thermodynamic limit we assume the standard 
fractal behavior, i.e. that the average mass belonging to 
a certain avalanche type inside a box of linear size £ is 
given by: 

M„(^;a) = M*(fTK'^° (33) 

in the limit ^ <C where ^ is the correlation length. The 
prefactor M*{a) is related to the concept of lacunarity 
|27l |. In general there can be fractals sharing the same 
fractal dimension, but with different lacunarities. The 
fractal dimension is related to the rate of change of the 
mass when the size of the box is changed. In contrast, 
the lacunarity is related to the size of the gaps of the frac- 
tal and is independent of the fractal dimension. In this 
way, the larger the typical size of the gaps, the higher 
the lacunarity. For many fractals, l28l| as lacunar- 
ity increases, the prefactor M*, decreases since the mass 
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FIG. 11: Collapses corresponding to h\{uV'^'^ ,vL^^^). 
The cuts of the scaling surface are taken at uL^^" = 
-0.58,0,1.22,2.7,3.8 and 5.0. Symbols correspond to the 
sizes indicated in the legend of the previous figure. The 
dashed line on the horizontal plane indicates the tendency 
of the crest of the scaling surface. The projection of the crest 
(Gaussian fit) is depicted with a continuous line on the plane 



inside a box of linear size € decreases. 

For finite systems it is necessary to translate the law 
into a finite-size scaling hypothesis. As usually done, 
we propose: 



e < L, (34) 



where the condition ^ ^ L stands for the fact that scaling 
only holds in the critical zone. Eq. H34() allows the data 
corresponding to the masses Ma to be collapsed and, in 
this way, to obtain the fractal dimensions da- We can 
predict the shape of Ma in two limiting cases: on the 
one hand, the scaling function Ma{uL^/^ ,1/ L) should 
behave as: 



L 



(35) 



in the limit ^/L<C^/L<lto recover the expression (|33f) 
from (P3)). On the other hand, in the Hmit £/L > ^/L, 
the scaling function corresponding to the subcritical 3D- 
spanning avalanche should behave as 



.iuL'^'',i/L)^ (\u\L'/^) 



I 



(36) 



if this avalanche fills a finite fraction \u 
for cr < (Tc in the thermodynamic limit 



of the system 
and we expect 



(37) 
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FIG. 12: (a) Number density of critical 3D-spanning 
avalanches at the critical amount of disorder obtained using 
Method 2. (b) Scaling plot corresponding to the data in (a) 
according to Eq. 11711 with 6 = 0.1. 



Such behavior can be obtained from Eqs. H3t)|) and H34|) 
using the hyperscaling relation /33_ = i/(3 — da-). 

Fig. El shows the scaled mass for all kinds of span- 
ning avalanches from simulations performed at a = cTc- 
We have only considered the mass of those SD-spanning 
avalanches that are equally classified by the two pro- 
posed methods. The fractal dimension rendering the 
best collapse is dj = 2.78 ± 0.05 for the ID-, 2D-, and 
critical 3D-spanning avalanches and ds- = 2.98 ± 0.02 
for the subcritical 3D-spanning avalanches. Both val- 
ues are in total agreement with those obtained indepen- 
dently in [23. Moreover, the slope of each collapse in 
the limit £ <^ L (left-hand side of the collapses) coin- 
cides with the fractal dimension used to obtain the col- 
lapses, so that the behavior (|35|l is confirmed. The pref- 
actors are M;{ac) = 0.95 ± 0.07, M^{ac) = 0.93 ± 0.07, 
M*^{ac) = 0.90 ± 0.07, and Ml_{ac) = 0.65 ± 0.07. The 
low value of the prefactor corresponding to the subcrit- 
ical 3D-spanning avalanches indicates that the gaps of 
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FIG. 13: (a) Number density of subcritical SD-spanning 
avalanches at the critical amount of disorder obtained using 
Method 2. (b) Scaling plot corresponding to the data in (a) 
according to Eq. 11711 with 9 = 0. 



these avalanches are large. As a consequence, the space 
filled locally by these avalanches is not as high as one 
would a priori think given the proximity of da- to 3. 

To study the behavior of M^- for a < cTc it is conve- 
nient to multiply the scaling function M^^ by the factor 
{e/L)~^. From Eqs. and it should behave as: 



/ (38) 

in such a way that, for a given value of uL^''^ < 0, 
the function {l/L)~^M^^ approaches a constant value 
for large values of £/L if the correlation length is fi- 
nite. Fig. [TKl shows the scaling collapses of {i/L)~'^M^- 
for three cuts of the scaling surface taken at uL^^^ ~ 
—3.73,-7.41 and —13.18. Note that such cuts are lim- 
ited from below at l/L = 1/80. In spite of this limita- 
tion, the results clearly indicate that (i) for small values 
of the behavior of {£/L)~^Mj,^ is power law with 
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FIG. 14: (a) Scaling collapses corresponding to the mass 
of the ID-, 2D-, and critical 3D-spanning avalanches at ct = 
(7c. The collapses corresponding to Mi and M2 have been 
displaced two decades upwards for clarity, (b) Scaling of M3- 
also at o" = o"c. In all cases the asymptotic behavior for small 
values of i/L has been indicated by continuous lines. The 
dashed line in (b) corresponds to the assimptotic behavior of 
Msc to compare it with Ma-. 



an exponent approaching da- — 3 ~ 0.02 (indicated by 
the dotted line) and (ii) for large values of £/L the func- 
tion {£ / L)~^ M3- tends to a constant value. (This latter 
tendency can only be observed for negative enough val- 
ues of uL^^") and confirms the hypothesis in Eq. 138|1 . 
Therefore, one can deduce that ^/L is finite for a < ac 
and it decreases when uL^/'^ becomes more negative. In 
addition, the results confirm the compact character of 
the subcritical 3D-spanning avalanche, as proposed in 
Ref. [23 by a different method. 



VII. DISCUSSION 

The results presented so far together with the results 
obtained in Ref. |23 provide a clear scenario for the phase 
diagram of the T = 3D-GRFIM with mctastable dy- 
namics in the thermodynamic limit. 
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FIG. 15: Collapses corresponding to 

{e.lLY^M-3,-{uL^''' ,ilL) for uL^I" = -3.73,-7.41 and 
— 13.18. The dashed line indicates the slope dz- — 3 = 0.02. 
Continuous lines are a guide to the eye. 



We have deduced that the subcritical BD-spanning 
avalanche occurring on tlie transition Une given by 
Eq. (|31|1 is compact and is thus responsible for the macro- 
scopic jump of the magnetization |2Cll |. Therefore we 
arc facing a standard first-order phase transition sce- 
nario with no divergence of the correlation length for 
o < Oc- At the critical point, this subcritical 3D- 
spanning avalanche becomes fractal at all length scales, 
and does not fill any finite fraction of the system. The 
end point (ctc, i^c) is a standard critical point. 

Fig. llGf a) shows the obtained phase diagram. The 
dashed line represents the first-order transition line given 
by Eq. (|31|l and the large dot the critical point. Note that 
this transition line is only approximate because it has 
been deduced from scaling arguments close to the critical 
point. Nevertheless it is remarkable that Eq. (|31|l for 
CT = renders (7J)3_(0) = —3.999 which is unbelievably 
close to the value —4 which can be computed by a (not 
so) trivial analysis of the coercive field of the hysteresis 
loop of the Gaussian 3d-RFIM with metastable dynamics 
corresponding to cr — *■ O+.fssjl 

In Fig. Iltif b) we also show for comparison the phase 
diagram of the 3D-GRFIM model in equilibrium at T = 
(ground-state) |23|. In addition, in Figs. El (c) and (d) 
we show the mean field (MF) solutions corresponding to 
both the metastable and equilibrium cases. 

The MF scenario indicates that the equilibrium and 
metastable critical points occur for the same value of 
CTc = z^^j-K (Table II V() . where z is the coordination 
number. In particular, 2 = 6 renders = 4.7873 
for the 3D-GRFIM. Below CTc, nevertheless, the transi- 
tion in equilibrium is a standard first-order transition 
(at H = 0), whereas it is a spinodal line Hgla) in the 
metastable case. From the equations in Refs. ^3 andfisl 



FIG. 16: Phase diagram corresponding to the exact 3D- 
GRFIM with metastable dynamics (a) and in equilibrium (b) . 
The phase diagrams corresponding to the mean field approx- 
imation of the 3D-GRFIM with metastable dynamics (c) and 
in equilibrium (d) are also shown. The thick continuous line 
in (c) indicates the spinodal transition in the metastable so- 
lution and the dashed lines in (a), (b), and (d) indicate the 
first-order transitions. In all cases we have indicated the co- 
ordinates of the critical point. 

it can be found that the metastability limit is: 



i?,(a)=a^21n^-z$err(^^ln^j, (39) 

when the external field is decreased. $err is the er- 
ror function. |3^ The continuous line in Fig. ^{c) cor- 
responds to Hg[a) for 2 = 6. The spinodal transition is 
characterized by a divergence of the fluctuations and the 
correlation length on the line Hs{(j) where the disconti- 
nuity of the order parameter Am occurs. 

In the present work we have shown that when com- 
paring the exact solutions (non MF) of both the equi- 
librium model and the metastable model, the character 
of the transition line docs not change. In both cases the 
transitions are standard first-order transitions with order 
parameter discontinuities and finite correlation length. 
This result agrees with the prediction that the tran- 
sition is abrupt for d < 8 as deduced from an e expansion 
analysis around d = 8. 

As indicated in Fig. [H and Table |IV| cr^'? = 2.270 ± 
0.004 in equifibrium and (7™*=* = 2.21 ± 0.02 in the 
metastable case. Regarding the value of the critical field, 
it is \Hc\ = 1.425 in the exact metastable model and zero 
in the exact equilibrium model. Thus, when the exact 
3d-GRFIM is studied, the critical point in equilibrium 
does not coincide with that corresponding to metastable 
dynamics. Nevertheless, the critical exponents are the 
same within statistical errors. The values are indicated 
in Table Hvl |3^ In fact, Dahmen et al. ^3 have already 
pointed out this similitude between the critical exponents 
for both models. These authors argue that agreement 
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between the two sets of exponents is rather unexpected 
since the two models are very different. Nevertheless, 
we can provide an argument based on rcnormalization 
group theory that indicates that the critical points in the 
two models (3D-RFIM in cquihbrium and the 3D-RFIM 
with metastable dynamics) correspond to the same fixed 
point in a more general parameter space. Within the 
framework of RG theory the critical surface (or critical 
line) is defined as the set of all points in the parameter 
space that flow to a certain critical fixed point when the 
renormalization group transformation is applied. The 
variation of the tunable parameters of a model describes 
a "physical" trajectory in the parameter space. Accord- 
ing to these definitions, the critical point corresponds to 
the point where the "physical" trajectory intersects the 
critical surface. The two models discussed here can be 
considered as particular cases of a more general model 
with the same 3D-GRFIM Hamiltonian and the follow- 
ing T = adiabatic dynamics: when B is varied, blocks 
of neighboring spins of size n < Umax flip when such a flip 
represents an energy decrease. The metastable dynamics 
introduced by Sethna corresponds to rimax = 1 (only sin- 
gle spin flips are considered) and the equilibrium model 
at T = (exact ground-state) corresponds to rimax = oo. 
The parameter rimax is a new parameter that must be 
considered in the RG equations. Since a critical point is 
found both with rimax = 1 and rimax = oo, it is natu- 
ral to assume that this is an irrelevant parameter. Thus, 
we propose the scenario presented in Fig. 1171 Changing 
n-rnax altcrs the position of the critical point, but not the 
critical exponents which correspond to the same critical 
fixed point. Numerical simulations of the 3D-GRFIM 
with rimax > 1 dynamics will help in clarifying this pic- 
ture. At present we guess that the RG flow follows the 
arrows schematically indicated in Fig. El Both the equi- 
librium critical point (ECP) and the metastable critical 
point (MCP) lie on the same critical surface (CS). In gen- 
eral a first-order phase transition occurs at the points in 
the parameter space that go towards a discontinuity fixed 
point when the RG transformation is applied. [33 We 
assume the existence of two discontinuity fixed points: 
the equilibrium discontinuity fixed point (EDFP) and 
the metastable discontinuity fixed point (MDFP). The 
EDFP controls the first-order phase transition in the 
equilibrium case {n^]^^ = and a < a^'^) and the MDFP 
controls the first-order phase transition when > 0. 

All the points that flow towards any of the discontinuity 
fixed points define the discontinuity surface (DS) where 
the first-order phase transition occurs. 

Another interesting question to be discussed is the de- 
termination of the correlation length ^ in the 3D-GRFIM 
with metastable dynamics. Avalanches can be under- 
stood as the zero temperature fluctuations in the driven 
system. Is their average linear size related to ^?. The 
flrst thing to note is that we have found that avalanches 
display two different fractal dimensions (and thus differ- 
ent associated /3 exponents). The v and fi exponents, 
nevertheless, are the same for all the scaling collapses. 



For instance, this is illustrated by Fig.|7|in this work and 
by Figs. 8, 9, and 10 in Ref. 0. Thus the behavior of 
the correlation length is unique: 

^^u-^Eiu^v-^'), (40) 

with S(a;) ^ a; -I- 0{x^). The fluctuations then can 
"choose" between two different mechanisms for propaga- 
tion either with fractal dimension df = 2.78 or with frac- 
tal dimension d^- = 2.98. The second point to be con- 
sidered is that ^ cannot be related to the size of the sub- 
critical 3D-spanning avalanche since we have found that 
^ is finite below CTc. Keeping these two observations in 
mind we propose that the correlation length is related to 
the average radius of the largest non-spanning avalanche. 
Below a = (7c the existence of a compact subcritical 3D- 
spanning avalanche does not allow for the non-spanning 
avalanches to overcome a certain finite length and thus 
^ is finite. Only at the critical point does the subcriti- 
cal 3D-spanning avalanche become fractal and allows for 
other spanning avalanches to exist and ^ becomes infinite. 
This behavior is much similar to what has been recently 
found in percolation. |^ [s^ We conjecture that some 
of the theorems that have been rigorously proven con- 
cerning the uniqueness of the infinite percolating cluster 
should be applicable to our case concerning the compact 
subcritical 3D-spanning avalanche. The present results 
should be considered as an interesting stimulus to pro- 
ceed with the analysis of percolation theory. For instance, 
we propose checking whether the fractal dimension of the 
spanning clusters is the same as that of the inflnite clus- 
ter at distances lower than the correlation length at the 
percolation threshold. 

VIII. SUMMARY AND CONCLUSIONS 

The results presented in this paper are mainly related 
to two topics in the 3D-GRFIM: firstly, the field de- 
pendence of the spanning avalanches, and secondly, the 
geometrical properties of the avalanches. We have ex- 
tended the FSS hypothesis proposed in to properly 
take into account the field dependence of the number 
densities na{H,a,L) and of the bivariate distributions 
Va{s, H]a, L). When carrying out such an extension, it 
is necessary to introduce a new scaling variable v and a 
new exponent fi related to the divergence of the correla- 
tion length when H approaches He {£, ~ {H — Hc)~^)- 
We have also introduced a scaling hypothesis for the field 
{H)a{<J, L) at which the different avalanches concentrate 
and their standard deviation (j^{(t,L). From the scal- 
ing collapses corresponding to the ID- and 2D-spanning 
avalanches we have found 1/n = 1.5. The study of 
the 3D-spanning avalanches is more intricate as already 
shown in , where we propose the existence of two dif- 
ferent kinds of 3D-spanning avalanches. In this paper we 
have proposed two approximate separation methods for 
classifying these avalanches as subcritical or critical. Us- 
ing these methods we have found that 1/n — 1.5 for both 
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Magnitude 


3D-GRFIM. Mean field 


3D-GRFIM 


Equilibrium Metastable 
(Ref. [29j) (Ref. [26]) 


Equilibrium Metastable 
(Ref. [22]) (This work) 


\Hc\ 


~ /o 1 ^ - /o 1 ^ 

ZyJ Z/TT Z\j Z/TT 




1.425 ±0.010 


1/m 


1/2 1/2 
3 3 
1/2 1/2 


1.37 ±0.09 1.2 ±0.1 
1.498 ±0.034 1.5 ±0.1 
0.017 ± 0.005 0.024 ±0.012 



TABLE IV: Coordinates of the critical point in the a — H plane and critical exponents for the 3D-GRFIM in equilibrium and 
with metastable dynamics corresponding to the MF approximation and the exact models. 

cr 



ECP 




FIG. 17: Schematic phase diagram of the proposed model in the space (_ff, n~ Jj^. , cr) in the adiabatic limit and T = 0. Arrows 
indicate the RG flow. The meaning of the different acronyms is explained in the text. The grey lines correspond to projections 
on the vertical planes. 



cases. Scaling enables the following behavior to be sck- 
etched in the thermodynamic limit: The ID-, 2D-, and 
critical SD-spanning avalanches only exist at the criti- 
cal point {(Tc, He) = (2.21, —1.425), where their number 
densities are infinite. In contrast, one subcritical SD- 
spanning avalanche exists below CTc and it occurs on the 
line (iJ)3-(cr) (Eq. (EB). 

From the average mass Ma{£; a) we have obtained the 
fractal dimensions corresponding to each of the types of 
spanning avalanches. This has allowed us to confirm in- 
dependently the results in df = 2.78 ± 0.05 for the 
ID-, 2D-, and subcritical 3D-spanning avalanches and 
(i3_ = 2.98 ±0.02. Furthermore, the behavior for a < (Tc 
of the mass corresponding to the subcritical 3D-spanning 
avalanches indicates that the correlation length ^ is finite 



below fTc. As a consequence, we conclude that the line 
(_ff)3_((T), where the discontinuity in the order param- 
eter occurs, corresponds to a standard first-order phase 
transition line and ^ only diverges at the critical point. 
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